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Abstract. 

We present 850pm imaging polarimetry of two high-mass star forming regions - W48 and S152. Within these 
regions we have identified two candidate high-mass protostellar objects - W48W and S152SE. The submillimetre 
continuum emission from the candidate HMPOs is bright in comparison to the nearby H II regions. W48W is 
a cold dense source, with no radio or mid-infrared emission. S152SE has an IRAS source IRAS 22566+5828 in 
the Southwestern part of the region, which appears in the mid-infrared 8.28/xm emission, but there is no radio 
emission. The 850/xm data shows another core within the region, in the Northeast. The polarimetry is ordered 
and the degree of polarisation is high over the candidate HMPOs - ~ 6% for W48W and ~ 8% for S152SE. 
Polarimetry results of this nature indicate a strong, ordered magnetic field threading the candidate HMPOs. The 
magnetic field direction in both S152SE and W48W is perpendicular to the direction of elongation of the cloud 
which would imply collapse along the field lines. Estimates of the magnetic field strength are derived using the 
Chandrasekhar & Fermi method. We calculate plane of the sky field strengths of ~0.7 mG for W48W and ~0.2 
mG for S152SE. We discuss the drawbacks of using the Chandrasekhar & Fermi method with a large beam size. 
Mass-to-flux ratios have been calculated and both clouds are found to be roughly critical. 
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1. Introduction 

Theoretical and observational studies of low mass 
star formation have lead to an evolutionary scenario 
based on the identification and studies of what are 

stellar objects (YSOs) 
Andre et all Il993l iBossI 
The precise details, such 
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earliest stages of star formation are of particular interest, 
as this is the stage where the highest rates of accretion 
occur along with the onset of outflows. The youngest 
protostars - the class YSOs - are somewhat difficult to 
observe since they are enshrouded with dust and so suffer 
large amounts of extinction. 

High mass star formation is not as clearly under- 
stood. Very few studies have been carried out on the ear- 
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liest stage s of high mass star fo r mation compared to low 
mass (e.g. iBeuther et alJl2002al lEvans et al]l2002jh This 
is generally because regions of high mass star formation 
are located at greater distances than the low mass star 
forming regions, and given that high mass star formation 
takes place on shorter evolutionary timescales, coupled 
with the statistical rarity of massive stars, it is difficult 
to locate and observe high mass stars in the process of 
formation - the so-called high-mass protostellar objects 
(HMPOs). Surveys searching for HMPOs are attemp ting 
to overcome th i s, in pa rticular lSridharan et alJ (j2002l) and 
IBeuther et al.1 l)2002bj) , who use a colour selecte d, radio- 
quiet IRAS sample, and iThomoson et al ] l|2002|) who are 
searching for HMPOs in the vicinity of ultra-compact H n 
(UCHu) regions. 

Mapping magnetic fields may yield vital information 
on the function of magnetic fields in star formation. 
Submillimetre polarimetry is one of the best ways to study 
the magnetic field. Spinning dust grains align themselves 
with the magnetic field, producing polarised thermal emis- 
sion which can then be measured and used to trace the 
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magnetic field structure projected onto the plane of the 
sky. 

In a bid to understand the magnetic field structures 
in the envelopes of YSOs we have obtained submillimetre 
imaging polarimetry of a sample of high and low mass 
YSOs using the JCMT in Hawaii. In this paper we present 
results from two high mass regions - W48 and S152. The 
two sources were selected as part of a larger sample of 
high-mass star forming regions. They were not selected 
specifically to find candidate HMPOs. 

2. Observations and Data Reduction 

The observations were carried out at the James 
Clerk Maxwell Telescope (JCMT) in Hawaii using the 
Submillim etre Common Use r Bolometer Array (SCUBA) 
camera (Hol land et al.lll999fl to gether with the polarime- 
try module ((Greaves et all2003h mounted on the entrance 
window of SCUBA. The data were obtained using the jig- 
gle mapping mode, and a 16 point jiggle map was com- 
pleted for each of the 16 different positions of the half- wave 
retarder (separated by 22.5°). This process was repeated 
8 times for W48, giving an integration time of 2.3 hours 
and 9 times (2.6 hours) for S152. The chop-throw was 150" 
with a position angle of 0° in local coordinates for W48, 
and 140" with a position angle of 90° in azimuthal coor- 
dinates for S152. The chop-throw size and position angle 
were chosen carefully to ensure that the telescope was not 
chopping onto other sources, contaminating the data. 

Data reduction was carri ed out using SCUBA User 
Reduction Facility (SURF) l|jenness fc Lightfootl fl99^) 
commands. The nod of the telescope and the flatfield cor- 
rections were carried out in the usual ways, with the ex- 
tinction correction usin g the Tsso calculated from the Tcso 
({Archibald et all l2002|) (measured to be ~0.08 through- 
out observations of W48 and ^0.07 for S152) for the time 
of the observation, calculated from polynomials fitted to 
measurements of tcso taken throughout each night of 
observing. Noise observations were used to identify and 
'switch off' any excessively noisy bolometers so that they 
took no further part in the data reduction. Sky removal 
of the data was carried out using bolometers with no sig- 
nificant flux from the source in order to remove sky noise 
from the data. 

Once the instrumental polarisation was calculated and 
removed on a bolometer-by-bolometer basis, the data were 
combined together to produce final maps in the three 
Stoke's parameters, I, Q and U using POLPACK software 
I Berrv fc C T ledhilll200lh . The Stoke's para meters were cal- 
culat ed by fitting the following curve ijSparks fc Axon! 
1999) to the data: 



I'k = -(/ + e(Qcos2<^. + [/sin20 fc )) 



(1) 



Where I' k is the expected intensity in image k, t is the 
wire-grid analyser transmission factor, e is the analyser 
polarising efficiency factor and 4>u is the effective retarder 



position angle after correction for the parallactic angle for 
image k. 

These were then combined to calculate the polarised 
intensity, P, and the position angles: 



p = VQ 2 + u 2 



n 1 U 

= — arctan — 
2 Q 



(2) 
(3) 



where is the position angle of the polarisation vectors. 
The uncertainties in each of these quantities are given 

by: 
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28.6° 



(4) 
(5) 



where s p is the signal-to-noise in P. 

The vectors were selected so that all have errors in po- 
larisation of less than 0.5% for W48 and and 0.75% for 
S152. The error limit was selected to be different for the 
two sources as the size of the SI 52 region meant that tele- 
scope offsets were required to cover the source, resulting in 
each pixel receiving less integration time relative to W48. 
The vectors are clipped on the polarisation errors instead 
of signal-to-noise as clipping on the latter would result in 
disposing of points where t he polarisation is low . 

Finally, flux calibration l|Jenness et alJl2002|) was car- 
ried out using observations of Uranus and CRL2688 for 
W48 and S152 respectively. The flux calibration factors 
(FCFs) on each night of observing were calculated to be 
2.02 Jy/arcsec 2 /V for W48 and 1.82 Jy/arcsec 2 /V for 
S152. These FCFs were then applied to the data. 

3. Results & Discussion 

3.1. Submillimetre Polarimetry 

FiguresHandHpresent the data observed with the JCMT. 
The greyscale is the Stoke's I, or the intensity of the ther- 
mal emission (850/im) from dust grains, the vectors repre- 
sent the direction of the magnetic field on the plane of the 
sky (the measured E vectors rotated through 90°). The 
length of the vectors indicate the degree of polarisation. 

The massive star forming region W48 (Fig. is 
a Galactic Hll reg i on at a distance of about 3.4 kpc 
llVallee fc MacLeod! Fmifflh and is ~ 1.3° to the East 
of the supernova remnant W44. Numerous infra-red 
sources have be en detected toward s W48, including 
IRS1 and IRS2 llZeilik fc Ladal Il978l). along with H?Q 
and OH masers llJenness. Scott fc Padmanlll995t rTurnerl 



1979HOenzel fc Downeslll977l: iHofner fc Churchwel]lll996l 
Caswellll200lF ~ 

Fig-fflshows that W48 has two submillimetre cores, the 
main, bright core in the East, W48main, and the smaller 
core in the West, W48W. The two cores appear to be con- 
nected by a ridge of dust and gas. In the Northwest of 
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Fig. 1. SCUBA polarimetry of W48. The greyscale is the 850 /im continuum emission, the overlaid vectors have been 
rotated through 90°to represent the direction of the magnetic field lines projected onto the plane of the sky. Epoch 
J2000. 



W48main the polarimetry B-vectors have a North-South 
direction, perpendicular to the ridge connecting the two 
cores. The degree of polarisation drops with increasing 
intensity, from ~2% towards the edge of the source, to 
~0.5% towards the intensity peak. T his has been observed 
in m any sources, i ncluding W51 fchrysostomou et all 
l2002h. Bok globules faenning et alJl200ll) and the OMC 
ijMatthews et al.ll2002[) . Explanations for this decrease in 
polarisation with increasing intensity include i) a twisting 
magnetic field towards the centre of the core, resulting in 
a lower net polarisation percentage - due to cancellation 
effects within the beam, ii) a lower grain alignment effi- 
ciency towards the centre of the core - producing more un- 
polarised emission, and iii) the dust particles may be more 
spherical in regions of high density, which also produces 
more unpolarised emission, resulting in a lower percentage 
polarisation. 



The polarimetry across the submm source W48W and 
the ridge connecting it to W48main is ordered (total vari- 
ance in the position angle, cr 2 ~ 0.05 radians 2 ) with a high 
(~6%) degree of polarisation. There is no evidence of a de- 
crease in polarisation with intensity towards the centre of 
W48W. 



Figure 121 shows the complex region Sharpless 152 
i SharplesslTl95flt) Loc ated at a distance of 5kpc 



i Wouterloot et al1ll993j) . S152 f orms part o f a physical 



pair with the nearby region S152 IjPismis 1990). Numerous 
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. The S152 Hn region is located 



in the Northwest of the imaged region in Fig. [2 where the 
ridge splits in two. The polarimetry in this region indi- 
cates that the magnetic field is North-South in direction 
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Fig. 2. Same as fig-E but for the region S152. 



(in the plane of the sky), aligning with the North-South 
ridge. 



In the Southeast of the imaged region there is a submil- 
limetre bright core S152SE. This is the first time that this 
source has been observed in the submillimetre continuum, 
with the S152 Hn region hav ing only been observed in 
the su bmillimetre previously bv lJenness. Scott fc Padmanl 
<ll995h . The polarimetry implies that the magnetic field 
lines are perpendicular to the direction of elongation of 
the source, suggestive of collapse along the field lines. The 
degree of polarisation is relatively high (^8%) across the 
region S152SE. In the Southwest part of S152SE there is 
an IRAS source (IRAS 22566+5828). 



3.2. Comparison with 1.4GHz and 8.28fim emission 

We have obtained 1.4GHz emission maps of the W48 
and S152 regions, downlo aded from the NRA O/VLA Sky 
Survey library (NVSS) ()Condon et al.lll99Sj) . Also, we 
obtained mid-infrared data of both of these regions at 
8.28^tm (A band) downloaded from the Midcourse Space 
Experiment (MSX) catalogue. Figs and 0] show contour 
plots of the 1.4GHz and MSX A-band emission overlaid 
onto our submillimetre measurements. 

The left panel of Fig. shows the W48 region im- 
aged with SCUBA, with contours of the 1.4GHz contin- 
uum emission. The radio emission aligns well with the H II 
region to the SE of the main submillimetre core. There 
is no radio emission associated with W48W. The right 
panel of Fig.|21shows MSX A-band (8.28 /im) mid-infrared 
data contoured over the SCUBA data. The infra-red emis- 
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Fig. 3. Same as figE but with contours of 1.4 GHz emission overlaid in the lefthand panel, a nd MSX A-band emission 
overlaid in the righthand panel. The stars represent the IRS sources iZeilik fc Ladalll978h. the diamond marks the 
position of the H n region (SIMBAD database), th e triangles represent water masers (Hof ner fc Churchwelfl 9961 , and 
the hexagon represents OH masers llCaswellll200l . The MSX data has a resolution of 20" , and the NVSS data has a 
resolution of 45" . 



sion occurs only over the main source, indicating that the 
Western source, W48W, is cold. 

Figure 0] shows the 1.4GHz emission (left panel) 
and MSX A-band data (right panel) contoured over the 
SCUBA data of S152. The 1.4GHz emission traces the Hn 
region, as marked by the white diamond. The MSX data 
reveal that there is infrared emission from the dust ridges 
that are close to the H n region, indicating that they are 
warm. There is also what appears to be a point source of 
warm emission in the submillimetre bright Southeastern 
source S152SE. This is probably associated with the IRAS 
source 22566+5828, marked by the white star as the offset 
could be due to the large error of the IRAS position, which 
encompasses the MSX source, or a temperature/optical 
depth effect. 

Comparisons of the SCUBA, MSX and NVSS 1.4GHz 
emission indicate that both W48W and S152SE are cold, 
dense regions of gas and dust. This suggests that they 
could be candidate HMPOs at an early stage of evolu- 
tion Most current HMPO surveys (e.g. ISridharan et alJ 
2002) use selection criteria which include only the more 
evolved HMPOs for study. W48W and S152SE are too 
cold to be included in these surveys, which are based on 
the IRAS Point Source Catalog, indicating that there may 
be many more young HMPOs which may only be inden- 
tified through submillimetre surveys. 

3.3. Masers 

Maser emission is common among star forming re- 
gions, specifically H2O, OH, methanol and ammonia 



maser emission. Maser emission is a good indicator 
of star formation as it only occurs in regions of high 
density. It has been shown that H2O maser emission 
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There are many masers within the regions observed. 
W48 has H2O and OH masers located around the 
W48main submm core. Whilst no maser emission has been 
detected on the source W48W, there has been a detection 
of an H2O maser close to the source, on the ridge con- 
necting the two cores. This indicates that the ridge and/or 
W48W is the site for early star formation, supporting the 
idea of W48W being a candidate HMPO. The location of 
the masers are indicated on Fig. 

Numerous masers have been detected towards S152, 
mainly H 2 and SiO, with H 2 masers having been stud- 
ied in S152SE. Again, only H2O masers have been discov- 
ered in this region indicating that the region is possibly 
undergoing the earliest stages of cloud collapse, again sup- 
porting the idea that S152SE is a candidate HMPO. Fig 
0] shows the location of the masers in this region. 

3.4. Masses of Candidate HMPOs 

The mass of the clumps can be calculated via the relation: 



M = 



gS v d 2 



K v B v (Tdust) 



(6) 
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Fig. 4. Same as figEO but the stars represent the IRAS sources 22566+5830 & 22566+ 5828, the diamond mark s the 
position of the Hn region (SIMBAD database), the triangles represe nt water masers J Compretto et al.||l99C). and 
the hexagons represent other masers including OH, methanol and SiO l|Szvmczak. et al . 2Q00tlHari u et alJll9 98). The 
MSX data has a resolution of 20" , and the NVSS data has a resolution of 45" . 



where g is the gas-to-dust ratio, S v is the flux density of 
the clump, d the distance to the clump, k v the absorption 
coefficient at frequency v, and B v (T,i ns t) is the Planck 
function at frequency v for a blackbody of temperature of 
Tdust ■ 

Using a gas-to-dust ratio of 100:1 (Hild ebrandlll983h . 
and an absorpt ion coefficient at 850 ixm of .15 m 2 kg -1 
estimated from lOssenkopf fc Hennina l)l994|) based on a 
number density of uh = 10 5 cm" 3 , thick ice mantles and a 
formation timescale of 10 5 years, temperatures of ~ 20K, 
flux densities and radii measured from the SCUBA data 
listed in tabled along with the distances quoted in table 
^ we calculate the masses to be ~ 1500 Mq and 3500 Mq 
respectively. 

The mass of these objects are subject to the un- 
certainties within the parameters used, specifically the 
gas-to-dust ratio, which could be as small as 45:1 
ijMcCutcheon et ail ^995). Also, the absorption coeffi- 
cient at 850 /zm has still not been determined precisely 
<|Hildebrandlll983l IChini et alJl200ll). The value adopted 
in this paper from lOssenkopf fc Henninel lll994h a grees 
with the values de termined by iBianchi et al.l l)2003^ and 
IVisser et all l|2002l) . 



3.5. Magnetic Field Strength 

It is important to try to establish the magnetic field 
strength as well as the magnetic field morphology in 
these regions for a more complete understanding of 
the contribution of magnetic fields to the star forma- 
tion process. Estimates of the magnetic field strength 
can be gained through use of the Zeeman effect and 
the IChandrasekhar fc Fermil l|l953[) (CF) method. The 
Zeeman effect uses line splitting to calculate the strength 
of the line-of-sight component of the field, while the CF 
method estimates the strength of the plane of sky field 
based on the dispersion of position angles of the polar- 
isation vectors. The CF method is applicable here and 
leads to estimates of 0.65 mG for W48W and 0.2 mG for 
S152SE. The parameters used are shown in tabled 

(B pos ) = fV^P— G (7) 

where p is the mean density (g cm -3 ), cr„ loa the line-of- 
sight velocity dispersion (cm s _1 ), erg is the dispersion 
in polarisation position angles and is corrected for mea- 
surement errors (of = o^ easured - (Terror) wher e is in 
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Object 


RA 


Dec 


Distance 


Flux 


Radius 


Mass 


Density 1 


°"«loi, 




ae 


-Bpos 




(J2000) 


(J2000) 


(kpc) 


(Jy) 


(cm) 


(M ) 


(g cm -3 ) 


(cm 


s" 1 ) 


(rad) 


(mG) 


W48W 


19 01 42.5 


+1 13 33 


3.4 


22 


9 x 10 1 ' 


1500 


1 x 10 _1B 


8 x 


10 4 


0.22 


0.65 


S152SE 


22 58 50 


+58 45 22 


5 


22 


2.3 x 10 18 


3500 


1 x 10~ 19 


8 x 


10 4 


0.26 


0.2 



Assuming spherical geometry. 



(rad), and / is a correction fac tor found to be 
(jo-striker. Stone fe Gamnfiell200lj) . 



0.5 



against collapse. Applying our magnetic field strengths to 
the Virial Theorem: 



The CF method should be used with caution when 
calculating the magnetic field strength with data of this 
resolution, as it can over-estimate the field strength. The 
large beam size of the JCMT means that small scale tan- 
gling of the magnetic field can occur within a beam, so the 
measured vectors only represen t the n et magnetic field 
within the beam. iHeitsch et alJ l|200ll) used their mod- 
els to test the effects of limited resolution on the ob- 
served field structure and the CF method of calculating 
the mean field strength. They conclude that small-scale 
variations were completely supressed, resulting in the im- 
pression of a more homogeneous field, and when this is 
the case, the magnetic field strength is over-estimated by 
a factor of 2 for field strengths consistent with molecular 
clouds. However, they purposefully leave out any infor- 
mation on the polarisation percentage p, which would de- 
crease where tangling of the magnetic field occurs. In order 
to avoid over-estimating the magnetic field, we have only 
used the candidate HMPO regions for our field strength 
estimates as the field is presumed a-priori to be homoge- 
neous, and there is no decrease in the degree of polarisa- 
tion which implies that on large-scales the magnetic field 
is not tangled. We have also appli ed the correction factor 
I Ostriker. Stone fc Gamrnie] 200lh of 0.5, consistent with 
the findings of IHeitsch et al I ll200lh . 



Uncertainties in our estimates of the magnetic field 
strength arise from calculating the density of the clumps, 
which incorporates the errors involved in calculating the 
mass. Also, errors in calculating the volume of the clump 
contribute as a spherical geometry has been assumed for 
each clump. The velocity of the gas within the clump in- 
troduces another error as a FWHM of ~ 2 kms -1 has 
been used b ut it may be anywhere between 1 kms -1 
and 3 kms -1 l|Brand et all200l[|Thompson fc Macdonaldl 
120031) . Measurement errors are also introduced by the an- 
gle although these are relatively small in comparison to 
the other errors stated. 



3.6. Cloud Support 

The virial theorem provides a gross method of analysing 
the stability of the clumps. We will simply estimate the 
values of the various energy terms to see how much sup- 
port is provided by bulk motions and magnetic pressure 



-3GM 2 



5R 



> 2 x 0.27MAw 2 + 0.1B 2 R 3 



(8) 



(Mc Kee et al.ll993ft where G is the gravitational constant, 
M is the mass of the clump, R the radius, Av the turbulent 
velocity of the gas (FWHM) and B the magnetic field 
strength in Gauss. 

Substituting values from table [T] into Eq.|3 it can be 
shown that the gravitational energy is of the order of the 
sum of the kinetic and magnetic energy terms, indicating 
the clouds are gravitationally bound or contracting. 

There are currently two models of star formation; star 
formation driven by ambipolar diffusion and that driven 
by turbulenc e. The ambipolar diffusi on model of star for- 
mation IjCiolek fc Mouschoviaslll994j) is based on subcrit- 
ical clouds that at first contract via ambipolar diffusion. 
When the central parts become supercritical, the magnetic 
field cannot support the core and it collapses dynamically. 
The turbulence model of star formation states that tur- 
bulence provides support against cloud collapse, as the 
magnetic field is too weak to provide support (the core is 
supercritical). 

To determine the importance of the magnetic field as 
a support mechanism, the mass to flux ratio can be used. 



(K) 

\ _ \ $ / actual 

x -~m — 

V $ / critical 

where we take the critical value to be: 



critical 



0.4 



(9) 



(10) 



llMestell p999). If the ratio is lower than the critical value 
(A < 1), the cloud is said to be subcritical, and the mag- 
netic field can support the cloud perpendicular to the field 
direction irrespective of external pressure. We calculate 
mass-to-flux ratios, in units of the critical ratio, of A ~ 
1.2 for W48W, using a magnetic field strength of 0.65 mG, 
and A ~ 1.3 for S152SE using a magnetic field strength of 
0.2 mG. This implies that both W48W and S152SE are 
close to critical. 

The value of the numerical coefficient in Eq. ^| is 
model-dependent, with Eq. 1101 representing the model 
for a spherical cloud inclusive of surface terms. iMestel 
(1999) also calculates a numerical coefficient of 0.2 for 
an extreme spheroidal model of eccentricity 1, whereas 
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iNakano fc Nakamural l)l978|) use a value of 0.4 for a disk- 
like model (ignoring surface terms), and ICrutcherl |2003) 
uses a value of 0.2 for a spherical model (ignoring sur- 
face terms). This leads to ranges in A of 1.2 < A < 2.3 
for W48W and 1.3 < A < 2.7 for S152SE. These calcula- 
tions of A are based on the plane of the sky magnetic field 
strength, B pos . If the magnetic field does not lie entirely 
in the plane of the sky in these regions, the actual mag- 
netic field strength will be greater than that calculated, 

and therefore the value of A will decrease. 

Previous studies of magnetic fields fsee ICrutcherll2003l 
for a review) have all indicated critical/supercritical star 
forming cores, with the only sub critical re s ults be ing for 
the diffuse ISM and HI clouds. ICrutcherl l|2003aj) stud- 
ied prestellar cores, and found that after corrections for 
geometry applied statistically, they were approximately 
critical. These are at roughly the same evolutionary stage 
as HMPOs, but are the low-mass equivalent. Our results 
are consistent with W48W and S152SE being very early 
stage HMPOs, and are found to be close to critical. 

4. Conclusions 

We present the first submillimctrc polarimetry data taken 
towards HMPOs, found in the vicinity of the W48 and 
S152 Hn regions. The polarimetry across W48 is compli- 
cated. There is a drop in the polarisation percentage with 
increasing intensity which may be explained by twisting 
of the magnetic field towards the centre of the core, or by 
relatively increased emission from unaligned dust grains. 
Polarimetry of S152 indicates that the plane of sky mag- 
netic field has a North-South direction across the Hn re- 
gion. 

Comparisons of our submillimctrc continuum data 
with 8.28^m and 1.4GHz emission indicate that W48W 
and S152SE are candidate HMPOs, at an early stage 
of evolution. The presence of water masers in these ar- 
eas supports this hypothesis. Both candidate HMPOs ex- 
hibit high degrees of polarisation (6-8%) and CF calcu- 
lations provide estimates of the plane of sky magnetic 
field strength of ~ 0.65mG for W48W and ~ 0.2mG for 
S152SE. 

The polarisation vectors are well aligned in both re- 
gions, indicating the magnetic field lines are perpendicu- 
lar to the direction of elongation and imply collapse along 
field lines. This polarimetry pattern has been predicted us- 
ing models of subcritical cores col lapsing under the e ffect 
of gravity and ambipolar drift fsee lPadoan et al.l200l|h al- 
though our estimates of A show that these clouds are close 
to critical and as such we are unable to rule out turbulence 
driven star formation models. 
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